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Abstract.— Native chickens play at least four significant roles for Thai people. In this study, the genetic 
variation of native chickens in northeastern Thailand was quantified based on the mitochondrial DNA 
HVS-I sequences. The phylogenetic relationship between these Thai native chickens as well as those of other 
geographic regions was also evaluated. A high level of genetic variation was found in the native chickens 
from northeastern Thailand, somewhat similar to that reported in Indian and Vietnamese populations but 
higher than African and Chinese native chickens. This likely implies a greater effective population size and 
longer demographic history and is consistent with the hypothesis that Thailand is within one of the places 
where chickens were first domesticated. Phylogenetic analyses revealed that Thai native chickens were 
divided into three lineages. The first lineage was closely related to G. gallus gallus, G. g. domesticus and G. 
g. murghi whilst the other two lineages were closely related to G. g. spadiceus. Thus, at least three 
subspecies, G. g. gallus, G. g. murghi and G. g. spadiceus, are likely to have been involved in chicken 
domestication in Thailand. 
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(Fumihito et al., 1994, 1996). However, 
recent molecular studies have suggested 
likely multiple origins of domestic chickens, 
and that at least three subspecies contributed 
to chicken domestication (Liu et al., 2006; 
Kanginakudru et al., 2008). 

In addition to searching for the origin of 
chicken domestication, quantifying the level 
of genetic variation is also important, since 
genetic diversity is critically important for 
species adaptability. For domestic animals, 
maintaining genetic diversity assures further 
breed improvement and selection. This is 
additionally important for globally 
cultivated/distributed species and those used 
for different functions (such as meat versus 
egg production), both of which are true for 
domestic chickens, where diverse adaptations 


INTRODUCTION 


Native chickens play four important roles 
for Thai people; namely they are used as 
food, in religious activity, as decoration, and 
for sport (Takashi and Kazunoba, 2010). 
Archeological evidence has indicated that 
chickens had already been domesticated in at 
least two sites, one in northern China by 
about 6,000 years BC and the other later in 
the Indus Valley by about 2,500 years BC 
(Tixier-Boichard et al, 2011). Early 
molecular studies originally indicated that 
domestic chickens were likely to have been 
derived from the continental red jungle fowl, 
Gallus gallus gallus, and that domestication 
took place in Thailand and/or adjacent areas 
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are required to local environmental 
conditions, pathogen resistance and 
development (e.g. resource allocation to egg 
or meat production), for example. Previous 
studies have shown that the genetic diversity 
of chickens was relatively high in areas 
where the origin of chicken domestication 
was speculated to have occurred 
(Kanginakudru et al., 2008; Berthouly- 
Salazar et al., 2010). 

Although Thailand has been speculated to 
be within one of the locations where chickens 
were first domesticated, only a small 
sampling of Thai chickens was included in 
the genetic analyses (Fumihito et al., 1996). 
A recent molecular study also ignored 
samples from the northeastern region of 
Thailand (Riztyan et al., 2012). In this study 
we, therefore, assessed the genetic variation 
of native chickens from northeastern 
Thailand and then examined the genetic 
relationship between Thai chickens and 
domestic chickens from other geographic 
regions based on DNA sequence differences 
in the mitochondrial DNA hypervariable 
sequence-I (HVS-I). 


MATERIALS AND METHODS 


Sampling, DNA extraction, Polymerase 
Chain Reaction (PCR) and Sequencing.— 
Domestic chickens were sampled from 21 
locations in 10 provinces in Northeastern 
Thailand (Table 1 and Fig. 1). DNA was 
extracted using a genomic DNA extraction 
kit (RBC BioScience, Taiwan) as per the 
manufacturer’s protocol and kept at —20 °C 
until use. The mitochondrial DNA HVS-I 
sequence was amplified using primers L16750 
(5-AGGACTACGGCTTGAAAAGC-3; 

Fumihito et al, 1994) and H522 (5- 
ATGTGCCTGACCGAGGAACCAG-3; Fu 


et al. 2001), using the PCR conditions of Liu 
et al. (2006). PCR products were checked by 
resolving over 1 % (w/v) agarose-TBE gel 
electrophoresis with ethidium bromide 
staining and uv-light transillumination, and 
cleaned using a HiYield™ Gel/PCR DNA 
Extraction Kit (RBC BioScience, Taiwan) as 
per the manufacturer’s protocol. Direct 
sequencing on both strands was performed 
commercially at Macrogen sequencing 
service (Seoul, Korea), using the same 
primers (L16750 and H522) as in the PCR 
but in separate reactions. 


Lao PDR 


Gulf of Thailand Ù 





FIGURE 1. Collection sites, within 10 provinces, for 
Thai native chickens in northeastern Thailand. Details 
of the sampling locations are given in Table 1. 
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TABLE 1. Sampling locations and number of samples of native chickens from northeastern Thailand use in 


the present study. 








Sampling location (Code) 


Sample 
size 


Nucleotide 
diversity (7) 


Number of 
haplotypes 


Haplotype 
diversity (h) 





Si sa ket Province (SR) 
Ban Donlee, Prangku District 
Ban Pimai, Prangku District 
All 
Maha Sarakham Province (MK) 
Ban Nong Wang, Mueang Maha Sarakham District 
Ban Nong No, Mueang Maha Sarakham District 
Ban Dindam, Mueang Maha Sarakham District 
All 


Kalasin Province (KS) 
Ban Huai Mek, Huai Mek District 
Ban Nong Chum Sang, Nong Kung Si District 
Ban Naku, Khao Wong District 





yap. 
n Nong Chiangcha, Phu Khiao District 


0.0147 + 0.0089 





7 1.000 


0.0191 + 0.0124 





5 1.000 


1.000 








t 0.0102 





Khon Kaen Province (KK) 
Ban Phai District 





Buriram Province (BR) 
Ban Na Pho, Phutthaisong District 
Ban Jaluay, Satuek District 
Al 


Nakhon Ratchasima Province (NR) 
Ban Dang Noi, Bua Lai District 
Al 

Roi-Et Province (RE) 

Ban Mek, Changhan District 
Ban Pho Thong, Selaphum District 








3 1.000 0.0213 + 0.0167 








Sakon Nakhon Province (SK) 
Ban Tad Phu Wong, Waritchaphum District 
Ban Kham Bo, Waritchaphum District 
Ban Nong Pan, Waritchaphum District 
Ban Non Na Kham, Waritchaphum District 
Ban Tao Ngoi, Tao Ngoi District 

















All 4 0.800 + 0.172 0.0117 + 0.0074 
Mukdahan Province (MH) 

Ban Na Ta Bang, Nong Sung District 11 

All 11 3 0.691 + 0.086 0.0072 + 0.0044 
Total 45 30 0.970 + 0.013 0.0156 + 0.0082 

Data analysis— The consensus DNA were estimated in Arlequin ver. 3.5 


sequences were aligned using CLUSTAL X 
(Thompson et al. 1997), with a final visual 
inspection. Representative haplotypes were 
deposited in GenBank under accession 
numbers KC817504 - KC817533. Haplotype 
diversity (fA) and nucleotide diversity (m) 


(Excoffier and Lischer, 2010). 

Phylogenetic relationships were analyzed 
based on the mtDNA HVS-I sequences of the 
30 haplotypes (from 45 samples) of Thai 
native chickens derived in this study plus 
those sequences previously published in 
Genbank (Table 2). Phylogenetic analysis 
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TABLE 2. List of sequences obtained from Genbank used in the phylogenetic analysis. 








GenBank accession 














Species/subspecies Location code Reference 
Gallus varius Bali, Indonesia D82912 Fumihito et al. (1996) 
G. gallus India GU902237 Singh and Kumar 
India GU902221 (unpublished data) 
Egypt EF586879 
Vietnam HM462085 Ramadan et al. (2011) 
Vietnam HM462217 Berthouly-Salazar et al. 
Sichuan Province, China AF512337 (2010) 
Sichuan Province, China AF512336 
Yunnan province, China AF512327 Liu, et al. (2006) 
Yunnan Province, China AF512077 
Sichuan Province, China AF512091 
Yunnan Province, China AF512108 
Guizhou Province, China AF512128 
Henan Province, China AF512139 
Yunnan Province, China AF512151 
TEE E E Hubei Province, China 7 AF512316 E E aaa 
G. gallus murghi India EU847798 Kanginakudru et al. (2008) 
India EU847791 
pE AA DIEN I EE AOAN e DS s AAEE AEE ADEA EE S India a EU847784 ERSA ASLA asso seine OSs PUNEA O EEn eee DES PE a AA D AAAA PES NEESS 
G. g. gallus Thailand D82900 Fumihito et al. (1996) 
Tokyo, Japan D82902 
Palembang, South Sumatra, Indonesia D82903 
Palembang, South Sumatra, Indonesia D82904 
Palembang, South Sumatra, Indonesia D82905 Berthouly-Salazar et al. 
Vietnam HM462088 (2010) 
Vietnam HM462090 
__ West Sumatra, Indonesia D82917 
Singaraja, Bali, Indonesia D82897 
West Java, Indonesia D82898 
Lanpung, East Sumatra, Indonesia D82899 
ae AB007718 Miyake (unpublished data) 
G. g. spadiceus Thailand D82907 Fumihito et al. (1996) 
Vientiane, Lao PDR AB007721 Nishibori et al. (2005) 
China: Yunan Province AF512154 Liu et al. (2006) 








was then performed using neighbor-joining 
(NJ), maximum parsimony (MP) and 
Bayesian inference (BI) analyses. The 
distance based NJ analysis was performed in 
PAUP* (Swofford, 2002) using the best-fit 
selected model. Bootstrap support was 
estimated using 1,000 replicates. MP 
analyses was also performed in PAUP*, 
using a heuristic search with 1,000 random 
addition sequence replicates, TBR branch 
swapping, and MulTrees effect. Bootstrap 
support was estimated for 1,000 replicates. 
BI analysis was performed using MRBAYES 


3.04b (Huelsenbeck and Ronquist, 2001). 
The best-fit model for Bayesian analysis was 
selected by hierarchical likelihood ratio tests 
implemented in MrModeltest (Nylander, 
2004). The general time-reversible (GTR) 
model (Rodriguez et al., 1990) with a gamma 
distribution shape parameter of 0.6624 and a 
proportion of invariable sites of 0.6729 was 
selected. The BI analysis was run for 
2,000,000 generations, with a sampling 
frequency of 100 generations. Tracer version 
1.5 (Rambaut and Drummond, 2004) was 
used for visual inspection of the point where 
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the log likelihood is stationary. Trees 
sampled before this point were discarded as 
burn-in. The remaining trees of two 
simultaneous runs were included in posterior 
probability calculations. 

In all cases the corresponding sequence 
from Gallus varius was used as the outgroup. 


RESULTS 


Mitochondrial DNA sequence variation.— 
A 545 bp fragment of the mtDNA HVS-I 
was obtained from 45 samples of domestic 
chickens from 21 locations in 10 provinces 
(populations) in northeastern Thailand. There 
were 34 base substitutions, comprised of 28 
transitions and six transversions. No 
insertion/deletions (indels) were detected in 
our sequences. A total of 30 haplotypes were 
identified. Seven of the nine populations that 
had more than two samples had a haplotype 
diversity of 1.00, with one population only 
having a single sample. Of the two other 
populations that from Mukdahan Province 
showed the lowest haplotype diversity 
(0.691) followed by that from Sakon Nakhon 
at 0.80. The overall haplotype diversity for 
the 10 populations (45 samples) was 
0.970+0.013 (Table 1). 

The nucleotide diversity in the 
populations differed between all nine 
populations that had more than two samples 
bar two (Chaiyaphum and Khon Kean) and 
ranged from 0.0072 (Mukdahan) to 0.0226 
(Roi Et) with an overall nucleotide diversity 
for the 10 populations (45 samples) of 
0.0156+0.0082 (Table 1). 








Phylogenetic relationship— All three 
phylogenetic analyses (NJ, MP and BI) 
revealed a similar tree topology with G. 
gallus forming two main clades (Fig. 2 — 4). 
Most of the samples belonged to clade I, with 


just G. g. bankiva from Indonesia forming 
clade II. Clade I was further divided into four 
subclades (I-1, I-2, I-3 and I-4) and the Thai 
G. gallus samples were found in subclades I- 
1, I-2 and I-3 but not in subclade I-4 that 
contained only G. g. gallus from Indonesia. 
Two subspecies of G. gallus, namely, G. g. 
domesticus from Indonesia and G. g.murghi 
from India also clustered with G. g. gallus in 
subclade I-1. G. g. spadiceus from China 
clustered with G. g. gallus in subclade I-2 
and G. g. spadiceus from Thailand and Lao 
PDR clustered with G. g. gallus in subclade 
I-3. The results strongly suggested that the 
taxonomic status of these subspecies needs 
further examination. 


DISCUSSION 


Genetic variation is important for species 
adaptability. For domestic animals, it is 
important for breed improvement and 
selection. Previous studies have revealed 
considerable levels of genetic variation in 
native chickens in regions where the origins 
of domestication were speculated to have 
occurred. High levels of genetic diversity 
(nucleotide diversity of 0.018) were found in 
the native chicken populations in India based 
on mitochondrial HVS-I sequences 
(Kanginakudru et al., 2008), and broadly 
similar results were reported from 
Vietnamese chickens with an observed 
nucleotide diversity of 0.013 (Berthouly- 
Salazar et al., 2010). This contrasts with the 
low genetic diversity reported from chickens 
in Africa (0.00745), which is consistent with 
the historical evidence that suggested the 
introduction of chickens into this region was 
far more recent at approximately 700 BC 
(Mwacharo et al., 2011). 
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Gallus varius 
G. gallus India 


G. gallus SR1 
G. gallus MK3 
G. g. gallus Japan 
G. gallus SRO 
G. gallus India 
G. g. gallus Thailand 
G. g. domesticus Indonesia 
G. gallus Hubei China 
‘ subclade I-1 
| G. g. murghi India 


G. gallus CP1 
OFG. galus MK4 
G. gallus BR1 
69g C. gallus CP2 
5 5G. gallus REI 
G. gallus SK3 
íG. gallus MH2, SR7, KS4 
. gallus MH1 
G. gallus MH3, KS1 
G. gallus MK2, KS2 
G. g. spadiceus Yunnan China 





Hi Sa |G. gallus Victnam 


G. gallus Sichuan China 


G. gallus Vietnam 


57 . gallus SR6 
P eG. gallus SR3 Clade I 
G. gallus BR2 
G. gallus MK1 
G. gallus KK2 subclade I-2 
64 G. gallus KS5 
G. gallus KS6, SK4 
5 IgG. gallus Sichuan China 
69 G. gallus Yunnan China 
"2G. gallus Guizhou China 
54 G. gallus Yunnan China 
98 66¢ G. gallus SR4 
G. gallus SR5 
58ļs4=G. gallus MK5, KK1, RE2 
53G. gallus SK1 
G. gallus SK2 
G. galus KS3 
G. gallus NR1 
G. gallus Yunnan China 
90 G. gallus BR3 ; 
98 4G. g. spadiceus Lao PDR 
64 G. g. spadiceus Thailand 
836G. gallus Yunnan China subclade I-3 
G. gallus Vietnam 
65 G. gallus India 
G. gallus Egypt 
57 G. gallus Sichuan China 
724G. gallus Henan China 





4 |a. g. gallus Indonesia | subclade I-4 


95 


i G. g. bankiva Indonesia Clade II 








0.01 


FIGURE 2. Neighbor-joining phylogenetic tree of chickens based on 545 bp of mitochondrial HVS-1 
sequences. Bootstrap values are shown above or near the branch. Native chickens from northeastern 
Thailand were labeled with a code representing the province where the samples were collected as SR, Si sa 
ket; SK, Sakon Nakhon; MK, Maha Sarakham, CP, Chaiyaphum; BR, Buriram; RE, Roi-Et, KK, Khon 
Kaen; KS, Kalasin; NR, Nakhon Ratchasima; MH, Mukdahan. The scale bar represents 0.01 substitution 
per nucleotide position. 
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Gallus varius 
G. gallus India 


. gallus SR1 
G. g. gallus Japan 
G. gallus SR2 
G. galuus MK3 
G. gallus India 
G. g. gallus Thailand 
G. g. domesticus Indonesia 
G. g. murghi India 
G. g. murghi India 
G. g. murghi India 
G. gallus China 
G. gallus MH2, SR7, KS4 
G. gallus MH1 
G. gallus RE1 
G. gallus MK4 
G. gallus CP1 
G. gallus BR1 
G. gallus CP2 
G. gallus SK3 
G. gallus Vietnam 
G. gallus Vietnam 
60 G. gallus China 
58 G. gallus SR3 
G. g. spadiceus China 
G. gallus Vietnam 
G. gallus Vietnam 
G. gallus MK2, KS2 
G. gallus MH3, KS1 Clade I 
G. gallus MK1 
G. gallus SR6 
s Geib stad 2 
G. gallus KS5 
G. gallus KS6, SK4 


subclade I-1 


G. gallus China 
G. gallus China 
G. gallus MK5, KK1, RE2 
G. gallus KS3 
87 58 |G. gallus NR1 
51 G. gallus SR4 

G. gallus SRS 

6G. gallus SK1 

G. gallus SK2 
G. gallus China 
G. gallus China 
G. gallus China 
89 G. gallus BR3 
94 4G. g. spadiceus Lao PDR 
G. g. spadiceus Thailand 
61)G. gallus China subclade I-3 
G. g. gallus Vietnam 
G. gallus Egypt 
54 G. gallus India 
G. gallus China 
G. gallus China 


86 


61 G. g. gallus Indonesia} sybclade I-4 





81 


54 


86 G. g. bankiva Indonesia Clade II 








FIGURE 3. Maximum parsimony phylogenetic tree of chickens based on 545 bp of mitochondrial HVS-1 
sequences. Bootstrap values are shown above or near the branch. Native chickens from northeastern Thailand 
were labeled with a code representing the province where the samples were collected as SR, Si sa ket; SK, 
Sakon Nakhon; MK, Maha Sarakham, CP, Chaiyaphum; BR, Buriram; RE, Roi-Et, KK, Khon Kaen; KS, 
Kalasin; NR, Nakhon Ratchasima; MH, Mukdahan. The scale bar represents one nucleotide change. 
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0.99 


0.88 
0.99 


0.01 


Gallus varius 
G. gallus India 
G. gallus SR1 
G. g. gallus Japan 
G. gallus SR2 
0.82 G. gallus MK3 
G. gallus India 
G. g. gallus Thailand 
0.9629 G. g. domesticus Indonesia 


D.9. G. g. murghi India subclade I-1 
G. gallus China 
G. gallus MH2, SR7, KS4 
A G. gallus MH1 
0.98 G. gallus RE1 
g— G. gallus MK4 
p.989. G. gallus CP1 
G. gallus BR1 
0.54 G. gallus CP2 
G. gallus SK3 


| G. gallus Vietnam 


G. gallus China 
0.99 G. gallus SR3 
0.87 G. g. spadiceus China 
——1).63 G. gallus Vietnam 
G. g. gallus Vietnam 
G. gallus MK2, KS2 
G. gallus SR6 
G. gallus MH3, KS1 
G. gallus BR2 
G. gallus MK1 
G. gallus KS5 
0.60 G. gallus KS6, SK4 
G. gallus KK2 


0.71 


0.99 0.63 
subclade I-2 


G. gallus China 
G. gallus China 
q. 53— G. gallus China 
0.69 G. gallus China 
G. gallus MKS, KK1, RE2 
G. gallus KS3 
0.97 [~ G. gallus NR1 
0.53} 0.95—G. gallus SR4 
G. gallus SRS 
0.50 G. gallus SK1 
G. gallus SK2 





G. gallus China 
D.56 G. gallus BR3 
1.00 1.00 G. g. spadiceus Lao PDR 
G. g. spadiceus Thailand 
0.97—-G. gallus China 
G. g. gallus Vietnam 
G. gallus Egypt 
0.91 G. gallus India 
G. galus China 
G. gallus China 


subclade 1-3 


0.94 Q g. gallus Indonesia| subclade 1-4 





1.00 G. g. bankiva Indonesia 








Clade I 


Clade II 


FIGURE 4. Bayesian inference phylogenetic tree of chickens based on 545 bp of mitochondrial HVS-1 
sequences. Bootstrap values are shown above or near the branch. Native chickens from northeastern Thailand 
were labeled with a code representing the province where the samples were collected as SR, Si sa ket; SK, 
Sakon Nakhon; MK, Maha Sarakham, CP, Chaiyaphum; BR, Buriram; RE, Roi-Et, KK, Khon Kaen; KS, 
Kalasin; NR, Nakhon Ratchasima; MH, Mukdahan. The scale bar represents 0.01 substitution per nucleotide 


position. 
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The genetic diversity of Thai native 
chicken has been reported previously, based 
on 98 autosomal SNP markers, to show a 
similar level of genetic variation to those 
from neighboring Myanmar (Riztyan et al., 
2012). From that reported SNP based 
analysis, both the Thai and Myanmar native 
chickens possessed a higher genetic variation 
than the commercial domestic chicken lines 
(Riztyan et al., 2012), as expected. Our 
results for native chickens from the northeast 
of Thailand, based on mitochondrial HVS-I 
sequences indicated comparable levels of 
genetic variation (average nucleotide 
diversity of 0.0156) to those reported from 
India and Vietnam. The results support the 
hypothesis that the area around and including 
Thailand was one potential location where 
chicken domestication occurred (Fumihito et 
al., 1996). 

Consistent with the high levels of 
observed mitochondrial DNA genetic 
variation, phylogenetic analyses revealed that 
the mtDNA haplotypes of village chickens in 
Northeastern Thailand were dispersed in 
several clades. Each population possessed 
genetically divergent mtDNA haplotypes. 
Similar results have been reported for 
Vietnamese chickens, where a high diversity 
within populations with diverse mtDNA 
haplotypes was found (Berthouly-Salazar et 
al., 2010). On the one hand, the results 
indicated a high rate of individual movement 
within the sampling areas. On the other hand, 
high levels of genetic variation in a single 
village (i.e. sampling location) is important 
for the population to maintain genetic 
variation to counter the effect of genetic drift 
(Berthouly-Salazar et al., 2010). 

Fumihito et al. (1996) speculated that 
Thailand and / or nearby areas were the 
location of the origin of chicken 
domestication. According to this study, the 
red jungle fowl, G. gallus gallus was thought 


to be the sole ancestor of domestic chickens. 
In contrast, Liu et al. (2006) demonstrated 
that domestic chicken could have originated 
from at least the three independent locations 
of (i) Yunnan, in South/Southwest China, (ii) 
the adjacent areas (Vietnam, Myanmar and 
Thailand), and (iii) India. This study also 
revealed that the same three subspecies of G. 
gallus (G. g. spadiceus, G. g. gallusi and G. 
g. murghi) were likely to have been included 
in the origin of modern domestic chicken, 
although when they were introgressed is not 
clear. Similar results obtained from Indian 
subcontinent chickens (Kanginakudru et al., 
2008). 

The results of the present study also 
support multiple origins of the current 
domestic Thai chickens. Phylogenetic 
analysis revealed that G. g. spadiceus and G. 
g. murghi were genetically closely related to 
Thai village chicken. On the other hand, G. g. 
bankiva from Indonesia was genetically 
distinct, suggesting that this subspecies may 
not have been involved in modern (maternal) 
chicken domestication. This has also been 
noticed in other studies (Fumihito et al., 
1996; Liu et al., 2006; Kanginakudru et al., 
2008). 

In conclusion, we found high levels of 
genetic variation in native chickens from 
northeastern Thailand, which could be the 
result of the high rate of chicken movement 
between villages. Phylogenetic analysis 
revealed that Thai native chickens were 
divided into three genetically divergence 
lineages that were associated with three 
different subspecies of G. gallus (G. g. 
spadiceus, G. g. gallus and G. g. murghi) and 
so suggests that at least these three 
subspecies were included in the maternal 
origins of the domestication of Thai native 
chickens, which has been observed 
previously in native chicken from India. 
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